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Abstract 

We look for differential equations of the form 



^ Cl {x)y {l) {x) = X n y(x), 



i=0 



where the coefficients {ci(x)}°°^ do not depend on n, for the generalized Jacobi polynomi- 
als {P%' ' M ' N (x)}™ =o found by T.H. Koornwinder in 1984 and for generalized Laguerre 
polynomials {L" ,M ' N (x)}°°_ which are orthogonal with respect to an inner product of 
Sobolev type. 

We introduce a method which makes use of computeralgebra packages like Maple and 
Mathematica and we will give some preliminary results. 
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1 Introduction 



In 1984 T.H. Koornwinder (see ||) introduced the polynomials { p&<P> M > N (x)} which are 

L J n=0 

orthogonal on the interval [—1,1] with respect to the weight function 
r(a + /3 + 2) M mWl , p 



2 «+/3+i r ( a + !) r ^ + i) 



1 + xf + M5(x + 1) + NS(x - 1). 



These generalized Jacobi polynomials generalize the Legendre type and the Jacobi type poly- 
nomials found by H.L. Krall. See for instance Q, [1C] and f7|. 

As a limit case Koornwinder also found generalized Laguerre polynomials \ L"' (x) > 

L J 71=0 

which are orthogonal on the interval [0, oo) with respect to the weight function 

' x a e~ x + M5(x). 



T(a + 1 

It is well-known that these generalized Jacobi and generalized Laguerre polynomials satisfy 
a linear second order differential equation with coefficients depending on n, but of bounded 
degree. See for instance ||. 

Now we are looking for differential equations of the form 

oo 

^2ci{x)y^\x) = X n y{x), 

where the coefficients {cj(x)}°^ are independent of the degree n. In |g, II and § several 
special cases were treated. Later more special cases were found by A.M. Krall and L.L. Little- 



john. See (Tj, § and [11]. 



These differential equations we are looking for are studied in spectral theory. A survey 
of orthogonal polynomials and spectral theory was given in jl]. A large number of references 
can be found there. 

In [|| J. Koekoek and R. Koekoek found a differential equation of the above type for the 
generalized Laguerre polynomials \L"' M {x) \ for all a > —1. This result generalizes the 

I J n=0 

results found for special cases of the Laguerre type polynomials. In fact we showed that this 
differential equation is of infinite order in general, but reduces to finite order for nonnegative 
integer values of a. See [Q] for more details and section 3 of this paper for more results 
concerning the coefficients of this differential equation. 

In H R. Koekoek and H.G. Meijer introduced the polynomials |L" ,A//,Ar (x)| which are 
orthogonal with respect to the following (Sobolev) inner product : 

oo 

<f,9> = r , ~~rr / x a e- x f(x)g(x)dx + Mf(0)g(0)+Nf'(0)g'(0), 
T(a + 1) J 

where a > —1, M > and N > 0. These polynomials { L"' M ,N (x)\ generalize the 

I J n=0 

polynomials \l% m (x)}°° since L% M '°{x) = L% M (x). 

^ J n=0 

In this paper we deal with the problem of finding a differential equation of the above 
type for both the generalized Jacobi polynomials \ P%>P' M,N (x) \ and the polynomials 

I J 71=0 

{L% M > N (x)Y 



I oo 
I n=0' 
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The differential equation found in || was computed by hand, without the help of comput- 
ers. This is nearly impossible for these other cases. We need computers to handle the very 
huge expressions we have to deal with. 



2 The method 

In this section we will describe the method we use to discover a linear differential equation 
satisfied by polynomials which are defined as a linear combination of classical orthogonal 
polynomials and their derivatives. These polynomials depend on x, the parameters of the 
classical orthogonal polynomials and on some extra parameters, say M and N. 

As an example we look at the polynomials \ L2J> ' fx) \ which are defined by 

I J n=0 

L% M ' N (x) = A L^(x)+A 1 -^L^(x)+A 2 ^L^\x), 

where the coefficients Aq, A\ and A 2 depend on n, a, M and N. Moreover, Aq, A\ and A 2 
are linear combinations of 1, M, N and MN. These generalized Laguerre polynomials are 
defined in such a way that L" ,0,0 (x) = lI^\x). See || and section 4 of this paper for more 
details. 

Now we start from the well-known (second order) differential equation for the classical 
orthogonal polynomials, which is 

xy"(x) + (q + 1 — x)y (x) + ny(x) = 

in the Laguerre case. If we substitute the new generalized orthogonal polynomials into the 
left-hand side of this classical differential equation we get an expression which can be seen as 
a polynomial in M and N. Now we add a linear combination of M, N and MN depending on 
y, y', y", etc. to the left-hand side of this classical differential equation and set this expression 
equal to zero. We hope that this is the new differential equation we are looking for. 
Our example leads to the differential equation 



MY j a l {x)y^{x)+NY j b l {x)y^{x) + 

i=0 i=0 

oo 

+ MN^2d(x)y ( ' i \x) +xy"{x) + (a + 1 -x)y'(x) + ny(x) = 0. 

Now we substitute the new orthogonal polynomials into this new differential equation 
and then the left-hand side becomes a polynomial in M and ./V with coefficients involving 
the classical orthogonal polynomials. Since this differential equation must be valid for all 
possible values of M and N, and M and N are supposed to be independent, all coefficients of 
this polynomial must be equal to zero. This gives us several equations in terms of unknown 
coefficients and of the well-known classical orthogonal polynomials. We remark that in case 
of dependence of M and ./V we have to deal with a completely different problem, since then 
we have a polynomial in one variable instead of two. We want the coefficients of y', y" , etc. to 
be independent of the degree n. So, if we substitute small values of n into these equations we 
get rather simple equations for the coefficients we are looking for. In general, these equations 
are huge expressions and are very difficult to handle. But here we use computers to do these 
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huge calculations. The mathematics used is very simple, but the formulas are too big to do 
this by hand. The calculations we did here were done by using Mathematica and sometimes 
Maple. Both computeralgebra packages can do these calculations easily, but they use a lot of 
computermemory. On a workstation with 8 megabytes of internal memory some calculations 
took several hours. 



3 The infinite order Laguerre differential equation 

In [Q] the following theorem was proved without the use of computers : 

Theorem 1. For M > the polynomials |-L"' M (x)} satisfy a unique differential 

I J 71=0 

equation of the form 



M ai(x)y^(x) + xy"(x) + (a + 1 - x)y'(x) + ny(x) = 0, 

8=0 



(1) 



where {cii(x)}^. are continuous functions on the real line and {cii(x)} c *L 1 are independent of 
n. 

Moreover, the functions {aj(x)}^ are polynomials given by 
a (x) 

(2) 

<■',(./•)--> (-i.r ""I . .11 . . \(a + 3)i-jX J , i = 1,2,3, .... 



n + a + 1 
n — 1 



For a / 0,1,2,... we have degree[aj(x)] = i, i = 1, 2, 3, . . .. This implies that if M > 
the differential equation (|l|) is of infinite order in that case. For nonnegative integer values of 
a we have 

degree[aj(x)] = i, i = 1, 2, 3, . . . , a + 2 

degree[aj(x)] = a + 2, i = a + 3, a + 4, a + 5, . . . , 2a + 4 

cii(x) = 0, i = 2a + 5, 2a + 6, 2a + 7, . . . . 

This implies that for nonnegative integer values of a and M > the differential equation (|l|) 
is of order 2a + 4. 

This differential equation was found by setting y(x) = L^ M (x) and by substituting 
small values of n in (|l|). Since the coefficients {ai(x)}°l 1 are independent of n this gives 
us a\(x), a,2(x), as(x), . . . explicitly. Then the general form of ai{x) was guessed and the re- 
sult was proved. Note that ao(x) does not depend on x, but does depend on n. From the 
three special cases (a = 0, a = 1 and a = 2) found by A.M. Krall and L.L. Littlejohn the 
general form of this coefficient could be guessed rather easily and be proved too. 

Later we discovered that the coefficients {ai(x)}°l 1 have the following interesting property. 
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Theorem 2. The coefficients {ai(x)}'^ 1 of the differential equation given by and (^) 
satisfy 

oo / 

sm7ra cc „ / 1 



, vr (a + 2)(a + 3) ^ a + 4 

For nonnegative integer values of a we have : 

E (l) < x ) = (-l) a+k a k (-x), k = 1,2, 3, 



x , a > —1. 



i=k 



(3) 



(4) 



Note that this theorem implies for nonnegative integer values of a : 

oo oo 

V^ai(x) = and ^iaj(x) = ( — l) Q+1 x. 
Proof. First we prove (0). Changing the order of summation we find 



E a *( a 



i=i 



« * (-1)^+1 / a + A / a + 2\ , ^ , 



Now we use the well-known summation formula 



2-^1 



a, 6 



T(c-a-b)T(c) 
r(c - a)r(c - b) 



c-a-b>0, 0,-1,-2,... 



to find 



^ (-1)' (a + 2\, 



l_g (-a-2)i(a + 3)i 



Hence 



oo 
i=l 



E 



(J + 

-a — 2, a + 3 

i + i 
r(j) 

T(i + a + 3)r(j - a - 2) 

« + A m 

{\j-ljT(j + a + 3)r(j-a-2y 



1 



, 3 = 1,2,3,.... 



i=o 



a + 1 



T(j + 1) 



j T{j + a + A)T{j-a-l)' 



(5) 
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Finally we use 



to obtain 



-a - 1) 3 



T(a + 4) (a + 4),r(-a - l)(-a - 1)/ 



x 



E 



-x 



r(a + 4)r(-a-l)^(« + 4) j 



1 



T(a + 4)r(-a - 1) I a + 4 



sm it z 



T(Z)T{1-Z) 7T 



r(a + 4)r(-a - 1) (a + 2)(a + 3) T(a + 2)r(-a - 1) 

1 sinvr(a + 2) _ 1 



sni7ra 



(a + 2)(a + 3) vr (a + 2)(a + 3) vr 



This proves (||). 

To prove (H) we take a € {0, 1, 2, . . .} and start with 



E j^y< x ) = EE Y^yr (J ' ^ I " ' " I (a + :5] ' 



i=kj= 



Now we use 




OO l oo oo 

E E = E E 

j=fc j=l j=l i=max(j,fc) 



to find 



OO I 



E(-t)j««M = E(-i) 



.7 + 



i=k 



3=1 



i(a + V 



x' 



E 



,. (i — k)\ \ i — j 

i=max(j,K) \ 



(a + 3) 



t-j. 
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For j = 1, 2, . . . , k we obtain 



_ ^ (~l) t+fc f « + 2 N 



(a + 3)j 



+k-j 



i=0 



r-, y- (-a - 2)j +fc _j(a + 3) i+k -j 
{ ' ^ i\{i + k-j)\ 



i=0 



and for j = k + 1, k + 2, . . . we find by using the summation formula (ph 



E 



(-1)* a + 2" 



(i -k)\\i- j 



(a + 3) 



E 

i=0 



(i + j — fc)! I i 



(a + 3)< 



■ (-a-2)^(a + 3)j 



(-1) 3 ' 
(-l) j 



2-^1 



-a — 2, a + 3 
j - k + 1 

r(j - k) 



Hence 



OO I 



E 



^(i-fc)l 



ai(x) 



r(j - fc + a + 3)r(j - fc - a - 2) ' 
(-a - 2) i+fc _j(a + 3)j +fc _ 



i=0 



E 

j=k+l 



+ fc - j)! 
r(i - fc) 



+ 



For the last sum we find 

OO 

E 

j=k+l 



J-l T(j -k + a + 3)T(j -k-a-2)' 



r(j - k) 



E 

3=0 



a+r 

J-l)" T(j - /c + a + 3)r(j - k - n - 2 ) 

j+fc+i r (i + 1) 



a + l 



= (-i) k x k+l Y, 

which equals zero since 



r(j + a + 4)r(j - a - 1) 

(-<* - r(j + 1) 

J= J 0" + *;)! T(j + a + A)T(j-a-l) 



(- x y 



i 



for j = 0,1, 2,..., a + l 



and 

So we have 



F(j - a - 1) 

-a — l)j + fc = for j = a + 2, a + 3, . . . and A; = 1,2, 3, 



OO -j 



+ l\ j ^ (-a - 2)t +fc _j(« + 3)i + fc-j 



i V 



i=0 



+ k-j)\ 
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The inner sum equals zero if k — j > a + 2. Hence 

k 



a 



+ l\ ,■(-«- 2) fc - J -(g + 3) 



j=max(l,fe— a— 2) 



(fc-j)! 
-a — 2 + A; — j, a + 3 + k — j 



(-1) 



fc+i 



E 



j'=max(l,A:— a— 2) 



a + l\ (a + 2 y 



k - j + 1 



(a + 3) k -;i{-xy x 



-a — 2 + A; — j, a + 3 + A; — j 
k-j + 1 



1 . 



Now we use the Vandermonde summation formula 

-n, b 



2F1 



1) = ~tt~" ) n = 0,l,2,... 

/ (c)n 



to find 



2-^1 



-a — 2 + A; — j, a + 3 + A; — j 
k-j + 1 



(—a — 2) a+2 -k+j 

(k-j+ l)a+2-k+j 



-1) 



a+2-k+j 



Hence 



00 -i 



j'=max(l,&— a— 2) 

x ^ fa + l\fa + 2 

k-j 



a + lWa + 2' 

i-ij U-j, 



(a + 3)*,-^ 



( - irl sG"-i)( 



! (a + fyk-jij 



(-l) a+k k\a k (-x). 



This proves (||). 



4 Some preliminary results for the Sobolev Laguerre polyno- 
mials 

In this section we look for a differential equation of the form 

00 00 



i=0 



i=0 



+ MN^a{x)y^(x) +xy"(x) + (a + 1 - x)y'(x) + ray (a;) = (6) 



for the polynomials 



y(x) = L% M ' N (x) = AoL£>(x) + M^\x) + A 2 ^L^(x), 
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where the coefficients Aq, A\ and A 2 are defined by 
A 



l + M( n + a ) + n( " + 2) - ( a + ^N h + a ) + 
\n-ll (a + l)(a + 3) In - 2 I 



MN 



n + a\ /n + a + 1 



(a + l)(a + 2) U-lM ra-2 



,,/n + a\ (n-1) /n + a\ 2MiV / n + a\ /n + a + 1 
Ai = M — — +iV| .1 + 



(7) 



n 



(a + 1) \n-l (a + l) : 



n 



A, 



N fn + a 



+ 



MAf /n + a\ /n + a + 1 



(a + l)ln-l/ (a + 1 



?i 



n — 1 



For details concerning these generalized Laguerre polynomials and their definition the reader 
is referred to || and ||]. 

Of course, since L° ,M '°(x) = L"> M (x), the coefficients {cii(x)}°l are given by (|2|). 

Although the general form is still an open problem so far, we know that the differential 
equation given by @ is not unique as in the case of the differential equation ([[]). We introduce 
the notation 

6q(0, a, x) = 



bo(n, a, x) = &o(l> a i x ) + /3o( ra ; oi, x), n = 1, 2, 3, . . . 
6j(a, x) = 6o(l, a, x)b*(a, x) + /3j(a, cc), i = 1, 2,3, . . . 
co(0, a, x) = 

co(n, a, ac) = bo(l, a, x) + 7o(n, a, 2), n = 1, 2, 3, . . . 



and 



t Ci(a,x) = b (l,a,x)c*(a,x) + ji(a,x), i = 1,2,3, 

Now we can prove the following theorem : 

Theorem 3. The polynomials { L° ,M ' N (x)\ satisfy the following infinite order differ- 

L J n=l 

ential equation : 

00 00 

J2b*(a,x)yV(x) + Mj2c*(a,x)y®( x )=0, (8) 



i=0 



i=0 



where 



and 



bt(a,x) = ~£(-lYr)(a + l) i - j x>, i = 0,1,2,... 



c*(a,x) 



-1)' 



■x l , i = 0,1,2,... . 



Proof. The proof is very easy and is based on the observation that 



00 / \ 



, n > 1, k = 0,1,2 



(9) 



9 



and 



n + a\ (— n)j 



n / (a + 1); 



To prove (g) we change the order of summation to obtain 

-l) j fi 



k = 0,1,2. 



(10) 



J2K(a,x)D i + k L^(x) = 



i=0 



i=0 j=0 

oo oo 



3=0 i=j 



EE 



(a + l)i^xW* +k L^(x) 
(a + l) i - j x>D i+k L£>(x) 

(a + l)iX j D i+ i +k L£Xx) 



f; ^±^* f; t±l X 3 D i+i+k L («) i3 



i=o "' j=o 

Now we use the definition of the classical Laguerre polynomial 



4 a) ( 



■r) 



n + a 
, n , 



iFi 



-n 



a + 1 



to obtain 

i=o 



n + oA ^ (— 1) 



n 

n + a 

v n , 
n + a 



ES^' E 



m—i—j—k 



j 

EE 



3=0 

oo oo 



, (a + l) m (m - i - j - k)\ 



3 =0 m=j 

oo 



m=i+j+k 

{-\y (-n) m+l+k x m 
jl (a + l) m +i+k (m - j)l 



E 



( — n )rn+i+k x 



m m 



m=0 (« + l) m+l +k ml f^ 
n + a\ (-n) i+k 



E(-!) J 



m 



\ n J(a + l)i +k 
Hence, by using the summation formula (g) we find 



J^b*(a,x)D i + k L^(x) 



n + a\j^(a + l)i (-n) i+k 



i=0 



n 



f Q i\ (a + l) i+k 



n + a\ (-n) fc 



"2-^1 



-n + k, a + 1 
a + fc + 1 



1 



y n J(a + l) k 

(n + a\ r(w)r(a + fe + 1) = (-ra) fc 

^ n ) (a + l) fe r(n + a + l)r(fc) nT(A;) ' 

which proves (|9]). The proof of ( JlOj ) is much shorter : 



E4(«,^ i+fc 4 a) (= 

i=0 



n + a 
n 



i=0 



m=i+k 



(a + l) m {m — i — k)\ 
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n + a 
n 



i=o l - 

oo m 



' (« + l) m +fc (m - i)! 



n + a\ (-1)' (-n) m+fc x m 

n )h*h i] (« + l)m+*(m-i)l 



n + a 



, n ) ^ (« + 1 )m+fc m! 
n + a\ (-n) fe 



i=0 



y n ){a + l) k 

By using these formulas (||) and (|i~0|) and by using the definition (0) of the coefficients Aq, 
A± and A2 we find 

00 00 00 

A Y / bU^x)D l L^\x) + A 1 Y / bU»,x)D l+1 Li a \x) + A 2 Y / bU^x)D l+2 Li a \x) + 

8=0 8=0 8=0 

00 00 

+ MA c*(a, x)D i L^\x) + MA X ^ c*(a, x)L> i+1 4°)(x) + 

00 

+ MA 2 4 (a, x)D l+2 L^ (x) 



i=0 



8=0 



n + a 



= -A x + (n - 1)A 2 + M I " ' " I, In 
This proves (||). 



8=0 



In the special cases a = 0, a = 1 and a = 2 differential equations of the form (|6|) are 
found too. In these three special cases of integer values of the parameter a we find a linear 
differential equation of formal order 4a + 10. By formal order we mean that for special cases 
(M = or N = 0) the true order might be lower. We give the results, but we will not give 
any proofs here. 

4.1 The special case a = 
If we take a = in (|2~|) we find 



ao(n, 0, x) = -n(n + 1), n = 0, 1, 2, . . . 



and 



oi(0, x) 
a 2 (0, x) 
a 3 (0,x) 
04(0, a;) 



3x — —x 

-2x + x 2 
1 , 



0, i = 5, 6, 7, . . . . 



11 



For the coefficients {/3j(x)}°^ we ™ ^^is case 



and 



P0(ll, U, X) 


— —n 2 (n 2 — Fl n — 1 9 3 

— ^2 v ~~ ' ' ' ' ' ' 


0i(O,x) = 





02(0, X) = 




03(0, x) = 


9 2 1 , 
-3-x + -x --x 


/9 4 (0, x) = 


2 + 3x - yl' + — X" 5 - — X* 


5 (O, x) = 


27 2 11 3 1 4 

-2x H x 2 x d + -x 4 

2 2 3 


6 (O, x) = 


r 2 9 3 1 4 
_ 5x 2 + _ x 3 _ _ X A 


7 (O,x) = 


4 3 1 4 

— x 3 + -x 4 
3 3 


08(0, X) = 


1 4 

"12 X 


0i(O,x) = 


0, i = 9, 10, 11,... . 



Finally, the coefficients {7«(x)}°^ turn out to be 



and 



7o (n, 0, x) = -^n(n 2 - l)(n + 2)(2n + 1), n = 1, 2, 3, . . . 



7i(0,x 
72(0, x 

73(0, x 

74(0, x 

75(0, x 

7 6 (0,x 

77(0, x 

7 8 (0,x 

79(0, x 

7io(0, x 
7i(0,x 



= 



~2 X 



5x z - -x A 

35 2 
2 



5 4 
24" 



28x z 



1 



14x J + -x x 



60 



56 



35 



1 



= -21x z + — x d - —x 4 + — x 



12 



12 



10 



Qx z - I2x 6 H x 4 x 



15 



1 



3x J x 4 + -x° 



12 



-x 



12 



-x 



60 

0, i = 11,12,13, ... . 



12 



Note that we have in this case : 

4 



10 



]T (h(0, x) = J2 A(o, x) = 7i(0, x) = o. 



i=l i=l i=l 

Hence, we have found the following linear differential equation of formal order 10 : 



^-MNx 5 y^°\x) + i-MiV(5x 4 - x 5 )y (9) (x) + 



+ 
+ 
+ 
+ 
+ 

+ 
+ 



y (8) (x) + 



-!^MiV(72x 3 - 45x 4 + 4x 5 ) - j^Nx 4 
\MN{2,Qx 2 - 72x 3 + 20x 4 - x 5 ) - ^iV(4x 3 - x 4 ) 
^MiV(-252x 2 + 224x 3 - 35x 4 + x 5 ) + ^N(-10x 2 + 9x 3 - x 4 ) 



V 



{7 \x) + 



y 



(6) (x) + 



— M7V(1680x 2 - 840x 3 + 75x 4 - x 5 ) + -N(-12x + 81x 2 - 33x 3 + 2x 4 )l y (5) (x) + 
60 6 



1 

24 



MiV(-420x 2 + 120x 3 - 5x 4 ) + 



+ j^N(24 + 36x - 150x 2 + 34x 3 - x 4 ) - ^Mx 2 



^MiV(15x 2 - 2x 3 ) + -N(-6 -2x + 9x 2 - x 3 ) + M(-2x + x 2 ) 
-^MNx 2 + ^N(2 - x 2 ) + ^M(6x - x 2 ) + x 



yW(x) + 
yW(x) + 
y"{x) + [1 - (M + l)x] y'(x) + 



+ 120 n 



MN(n 2 - l)(n + 2)(2n + 1) + 10iVn(n 2 - 1) + 60M(n + 1) + 120 y(x) = 0, 



for the polynomials 



y(x) := I% M > N {x) = A L n (x) + A 1 L' n (x) + A 2 L" n {x) 



with 



and 



L n (x) :=4°)(x) = X: 



k=0 



(-l) k /V 



x k , n = 0,1,2,... 



A = 1 + Mn + -7Vn(n - l)(2n - 1) + — MNn 2 (n 2 - 1) 
6 12 

Ai = M + 7Vn(n - 1) + -MNn(n 2 - 1) 

3 

A 2 = Nn + ^MNn(n + 1). 



4.2 The special case a = 1 

In the special case that a = 1 we find successively 

1 



a (n, 1, x) = -n(n + l)(n + 2), n = 0, 1, 2, . . . , 



13 



and 



(3 (n, l,x) 

P 1 {l,x 
f3 2 {l,x 

(3 5 (l,x 
/? 6 (l,x 
(3 7 (l,x 
Ps{l,x 
(3g(l,x 



(3i(l,x 



ai(l, x 
0,2(1, x 

a 3 (l,x 
04(1, x 
a 5 (l,x 
a 6 (l,x 

CLi(l,X 
1 



—x 

6x — x 2 

-lOx + Ax 2 - -x 3 
6 

5x — 5x 2 + -x 3 

2x 2 - \x 3 

-x 3 
6 



240 

= 
3 
2 



= 0, i = 7,8,9, ... , 
n(n z - l)(n + 2)(3n - 1), n = 1, 2, 3, . . . , 



rX 



= -8 x + Ax z x 6 

2 12 

25 77 2 47 o 7 4 
= Y +8X -^ X + T2 X -48 X 



25 



77 



51 



23 



-6 x H x —x H -x x 



227 

6x - 34x z + —x 3 



, 2 

79 
IT 

79 4 1 5 

48 X+ 8 X 
1 , 



I l - —x* + — x 4 
6 6 



r 



24 

19 4 1 5 

¥ + 16 



-x 



80 



' 4 
-x X 

] 16 



80 



7o(n, l,x) 

7i(l; x 
72(l,a; 

73(l,a; 

74(1,2; 

75(1,2; 



= 0,i= 11,12, 13, ... , 

n(n 2 - l)(n + 2)(n + 3)(5n 2 + lOn + 2) 
10080 : 



n = 1,2,3,... 



"4 X 



5x 2 



115 



-35x z + — x" 



23 

48 : 



2 105 3 65 4 

119x 2 x 3 H x 4 - 

2 12 



31 

240' 



-x 



14 



76(1, X 
77(1, X 

79(1,2; 
7io(l, x 

711(1,2; 
712(1,2; 
713(1,2; 

714(1, X 



441 



49 



= — — x z + 147x d - — x 4 + 



-x 



228x z - 232x J + — x 



117 4 14^ 5 



29 

24 : 

■x L 



72' 



-x 



495 



837 



4 4 



1289 4 467 = 
"IT* + 48 * 



2016 

3 6 1 7 
-x D H x' 

8 288 



55 2 605 o 129 4 571 = 25 fi 1 - 

— x z x 6 H x 4 —x s H x b x' 

2 6 2 48 36 96 



121 o 671 4 257 = 55 

x x H x x u H x' 

6 24 30 72 288 



.6 



61 „4 27 . 



1 



6 



— x x H — x 

12 8 2 

9 5 13 6 1 - 

—x 5 - — x b H x' 

16 72 96 



288' 



-x 



36 



-x 



288 



2016 

0, i = 15,16,17,... . 



Hence, in this case we have 



10 



14 



£ Oi(l, x) = £ a;) = E x) = 0. 



i=i 



i=i 



i=i 



This implies that we have found a linear differential equation of formal order 14 for the 
polynomials I L^ M ' N (x) \ 

4.3 The special case a = 2 

If we take a = 2 we find successively 

a (n, 2, x) = ^n(n + l)(n + 2)(n + 3), n = 0, 1, 2, . . . , 



ai(2, x 
a 2 (2,x 

a 3 (2,x 

04(2, x 

a 5 (2,x 

a 6 (2,x 

07(2, x 

a 8 (2,x 
Oj(2,x 



= —x 



3 n 

lOx - -x 2 
2 

-30x + 10x 2 - ^x 3 

35x _ l^ x 2 5 3 _ 1 4 

2 2 24 



9 



1 



= -14x + 21x z x 6 + -x 

2 6 

7 1 

= -7x 2 + -x 3 x 4 

2 4 

= -x 3 + ix 4 
6 



24 

0, % = 9, 10, 11,... , 



15 



Po(n, 2, x) = j^n(n 2 - l)(n + 2)(n + 3)(2n - 1), n = 1, 2, 3, . . . , 



7o (n, 2,x) = 



/3i(2,x) 
/3 2 (2,x) 

/3 3 (2,x) 

/? 4 (2,x) 

/3 5 (2,x) 

/3 6 (2,x) 

/3 7 (2,x) 

/3 8 (2,x) 

/3 9 (2,x) 

/?io(2,x) 

/3n(2,x) 

/5i 2 (2,x) 
ft(2,x) 

ra(ra 2 




2 



6' 



-x 



50 25 2 7 , 

2x H x 2 x 6 

3 6 18 



50 



92 



11 



45 + —x - — x z H x 3 -x 4 

3 3 2 24 

290 n 

-49 - 45x + — ^ 

56 , 887 2 188 , 91 4 97 . 1 fi 

h 49x x 2 H x d x 4 H x 5 — x° 

3 6 3 12 360 540 



„ s 49 
27x 3 + — x 



13 
360' 



-x 



56 



217 



451 



169 



29 



-x 



2 

2 271 , 

Y x + IT 



92 



97 



173 



3 4 173 5 
2* + 360* ' 



45* +9Q 1 



540 

0, i = 13,14,15,. 



x H x' 

6 12 

337 4 5 5 1 6 

x 4 + -x 5 x 6 

24 4 36 

77 5 1 6 

—x 5 + —x 6 



1 



36 X + 90 X 



72 



27 



36 



x 



l)(n + 2)(n + 3)(n + 4)(2n + 3)(7n 2 + 21n + 2) 
1088640 ! 



n 



1,2,3, 



= 



— x 
6 

35 o 13 q 

—x 2 x 3 

6 18 

_ 70x 2 + _ 7 4 



413x^ 



49 



13 



161x H x x 



-1386x 2 + 



2317 



30 

245 4 119 . 73 f 

x H x x 

2 18 720 



21 



= 2827x 2 - 2189x 3 + 497x 4 - 42x 5 + — x 6 - 



16 



59 
5040 



x 



-3575x 2 + 3872x 3 



9779 4 , 29 fi 19 7 11 c 

x 4 + 148x 5 - —x 6 + —x 7 - — x s 

8 4 144 17280 



16445 2 77 935 , 11473 4 5797 . 
-x 2 rr— x 3 H — x 4 —x b + 



6 



18 



18 



3257 
144 ; 



31 

48 



x 



16 



11 8 1 Q 
J x X 

1728 77760 

, , 7007 2 17875 » 23023 4 4521 = 2651 fi 1303 7 

710 (2, x) = x 2 H x 6 x 4 H x & x b + — — x' + 

' uv ' ; 6 6 12 10 60 720 

121 o In 

x s -I x 

4320 8640 

, , 637 2 3458 , 7189 4 1235 = 13511 fi 2311 7 

7n (2, x) = x l x H x 4 x H — x b - — — x r + 

' v ; 3 3 6 3 240 720 

77 r 1 n 

J x X 

1080 2160 

, x 1729 , 1274 4 2119 . 845 fi 899 7 1001 „ 7 q 

7i 2 (2, x) = x 6 x 4 H x 5 x b H x 7 x 8 H x 9 

' V ; 9 3 9 18 240 8640 6480 

,0 \ 196 4 w = 1189 B 415 7 539 o 7 

/ x -5 361 fi I 423 7 77 o 7 

714 2, x = llx 5 - —x 6 + x 7 -x 8 + — — 

' v ' 1 48 1008 864 4320 



-x 9 



,x 9 



/ \ 361 a 2 7 11 o 7 n 

7i 5 (2, x) = —^--^—a*-—^ 

7l6(2 ' x) = Yo x7 ~ Tt1o x8 + 21V 



11 8 1 
r S 



X° - —TTX- 



717(2 ' X) = 8oio 
718(2, x) = ^x 9 
7i (2,x) = 0, i = 19,20,21,.... 



Hence, we have 

8 12 18 

]T <n(2, x)=J2 A(2, x) = J2 7i(2, x) = 0. 

2=1 2 = 1 2=1 



t J_ t J_ t _L 

This implies that the polynomials \l^ m,n (x)\ satisfy a linear differential equation of 

I J n=0 

formal order 18. 

5 The generalized Jacobi polynomials 

In this section we will deal with the problem of finding a differential equation for the gener- 
alized Jacobi polynomials \ P"'P ,M ' N {x) \ 

I J n=0 

Since the well-known second order differential equation for the classical Jacobi polynomials 

fp( a 'Ph x \\ i s given by 

I J n=0 

(1 - x 2 )y"{x) + \p-a-(a + P + 2)x] y'(x) + n(n + a + (3 + l)y(x) = 0, 
it is clear that we look for a differential equation of the form 

CO 00 00 

M (x)y^ (x)+NJ2 k (x)y^ (x) + MN £ * {x)v {i) (x) + 
i=0 i=0 i=0 

+ (1 - x 2 )/(x) + \p-a-(a + P+ 2)x] y'(x) + n(n + a + (3 + l)y(x) = 0, 



17 



where 



and 



y(x) = P%>P> M > N (x) = A Pt P \x) + [^(1 -x)- A 2 (l + x)] ^P^ip) 



A Q = 1 + M- 



n + (3\(n + a + (3 + l 

In + a\ 
\ n I 



+ N- 



n + a\ I n + a + (3 + 1 

K n - 1 )\ 



+ 



+ MN 



(a + /3 + 2) 2 /n + a + /3+l N 



(a + !)(/? + 1)1 n-1 



M 



n + /3\/re + a + /5 N 
nil n 



(a + /?+ 1) 



A, 



/n + a \ 

V n ) 

n + a\ (n + a + (3 
n \ n 



+ 



MN fn + a + p\(n + a + P + l 



a + 1)\ n — 1 



n 



(a + /3+ 1) 



n + /f 



+ 



MN n + a + 0\ n + a + (3 + l 



09+ 1) V 



71 



and 



Here we used the same definition as in Q, but in a slightly different notation. 

In this case the differential equation will be unique in its general form if it exists 
introduce the notation 

ao(x) = ao(n,a, (3,x) 

Oi(x) = ai(a, f3,x), i = 1, 2, 3, ... , 
b (x) = b (n,a,/3,x) 

bi(x) = bi(a,(3,x), i = 1,2,3, ... 

c (x) = c (n,a,(3,x) 

Ci(x) = Ci(a,f3,x), i = 1,2,3, 

Since the polynomials \ P?'P> M ' N (x)\ satisfy the symmetry relation 

I J n=0 

pa,f>,M,N(_ x) = ^npP,a,N,M^ 

we have 

ao(n, a, /?, — x) = &o(n, /3, a, x) 
di(a, (3, -x) = (-l) l bi(P, a, x), i = 1, 2, 3, . . . 
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and 



co(n, a, f3, —x) = co(n, 0, a, x) 



Ci(a, (3, -x) = {-l) % Ci([5, a,x), i = 1, 2, 3, 

The general form is still an open problem, but the symmetric case a = (3 and M = N 
seems to be much less difficult. All special cases known so far (see for instance [jl|) point out 

that if a = (5 and M = N we can choose = for all i = 0, 1, 2, Therefore, we will 

consider this special case first. 

5.1 The symmetric generalized Jacobi polynomials 

All known examples seem to point out that there might be a differential equation of the form 

oo 

M^ai(x)y (i) (x) + (1 - x 2 )y"{x) - 2(a + l)xy'(x) + n(n + 2a + l)y(x) = 0, (11) 

i=0 



where 



and 



2/0*0 = i£ 



(12) 



C = l + M^ n + 2a + 1 n + 2a + 1 

(a + 1) \ n / \ n - 1 



(13) 



Ci 



2M /n + 2a\ 2M 2 / n + 2a \ / n + 2a + 1 



(2a + 1)1 n J (a + l)\n-l ) \ n 
This differential equation turns out not to be unique. We write 

ao(x) := ao(n, a, x), n = 0, 1, 2, . . . 
aj(x) := ai(a, x), i = 1, 2, 3, 



If we substitute (|12| ) in the differential equation (11) then we finally find three equations 



for the coefficients {oi(x)}^ which are equivalent to the following two : 



i=0 



(2a + 1) 



n + 2a\ d 2 
n } dx 2 



Pt a \x) 



and 



We introduce the notation 

ao(n, a, x) = ao(l, a, z)6o(«j a, z) + co(n, a, x), n = 0, 1, 2, . . . 

cij(a, z) = ao(l, a, x)bi(a, x) + Cj(a, x), i = 1, 2, 3, . . . . 
Now we can prove the following theorem. 
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Theorem 4. The polynomials { p^ a > M M ( x \ X satisfy a linear infinite order differential 

I J n=0 

equation of the form : 

oo 

5>(x)y«(x) = 0, 



i=0 



where 



1 



bo(x) := bo(n, a, x) = - [1 - (-1)"] , n = 0, 1, 2, 



ji-i 



6,-(x) := 6,;(a,x) 



-(-x) 1 , i = 1,2,3, 



Proof. To prove this theorem we have to show that 

oo 

J2h(x)D l Pt' a \x) =0 



i=0 



Y j ib i {x)D i P^ a \x)+xY j b i {x)D i+1 P^ a \x) = 0. 



i=0 



i=0 



To do this we note that the polynomials { Pi a ' a \x)\ are defined by 

I J n=0 



n + a 



n 



I 2^1 



-n, n + 2a + 1 
a + 1 



1 — x 



Further we write 

OO -I oo 

^(x^'P^x) = - [1 - (-l)»]P^)(x) +^6 l (x) J D^)(x) 

i=0 i=l 

and 

OO -I j oo 

^(x^^P^'^x) = - [1 - (-1)"] ^Pt' a) (x) +J2bi(x)D i+1 Pt ,a) (x)- 



1=1 



Now we easily obtain 

oo 

][>(x)I^ Q )0z 

i=l 



1 ^n + a\ ~ x* ^ (_ n ) fe ( n + 2a + l) fc ^1 -x^ fe 



2 V " ;S »! ^ (a + l) fc (fe-i)! V 2 
1 /n + a\ ~ (-n)fe(ra + 2a + l) fc * / k\ , t fl-x\ k - 1 



i=l 



1 + x\ fe fl-x 



1 (n + a\ ^ (-n) fc (n + 2a + l) fc 
H n J fe tt (a + l) fc fc! 

Pi a ' Q) (-x)-P^ Q )(x)] =~ [!-(-!)"] P^(x), 
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1 I n + a 



i=0 



n 



E 



E 



-n)k{n + 2a + l) k / 1 — x 



k—i 



1 / n + a\ ^ (-n) fc (n + 2a + 1)^ ^ (k - l\ , 
n (« + !)*(*-!)! U ' '' 



2 v 
1 / n + a 



n 



E 



i - l / 

{-n) k (n + 2a + l) k (l + x^ k " { 
\ (a + l) fc (fc-l)! 



1 — x 



d 



{-l) n x^-P^ a \x) 



and 



Y,bi{x)D i+1 P^\x) 
i=i 

1 / n + a 



n 



OO j oo 
j=l fe=i 



(-n) fe+ i(n + 2a + /l - x 



(a + l) fc+ i(fc-i)! 



fc— i 



1 / n + a 



n 



1 / n + a 



n 



1 / n + a 



n 



1 (n + a 



E 

k=l 

oo 

E 

fc=i 

oo 

E 

k=0 

oo 



(-ra) fc+ i(ra + 2a + 



E 



E 



(a + l) fc+ ifc! , . 

(-re) fc+ i(n + 2a + l) fc+ i 
(a + l)fc+ife! 

(-n)fc + i(n + 2a + l)fc+i 
(a + l) fc+ ifc! 

(-n) fc (n + 2a + l) fc 



1 — x 



1 + x 
2 

1 + x x fc 



1 — X 
2 

1 x N 



- (a + l) fc (A;-l)! 



l + xX^ 1 (\ 



x 



fc— l 



= -I [1 + (_!)«] A P (^) (x) . 

This proves the theorem. 

For the coefficients {cj(x)}°^ we find by using packages like Maple or Mathematica 



co(0, a, x 
co(l, a, x 
co(2, a, x 
co(3, a, x 
Co (4, a, x 

cq(5, a, x 






4(2a + 3) 

4(2a + 3)(2a + 5) 

2(2a + 3)(2a + 5)(2a + 6) 

-(2a + 3)(2a + 5)(2a + 6)(2a + 7). 



Hence, we might guess that 

2n(n-l) /n + 2a + 2 



c (n,a,x) 
If we write 



n 



)=4( 2a + 3 )(" + 2 _° 2 +2 ),„ = 0,l,2„... 



(a + 2) V 

Cj(a, x) = (2a + 3)(1 — x 2 )c*(a, x), i = 1, 2, 3, . 
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we obtain 



c\(a, x) 
cZ>(a,x) 

c\(a,x) 
c%(a,x) 
c* 5 (a,x) 
c* 6 (a,x) 
Cj(a, x) 
c* 8 (a,x) 

c* 9 {a,x) 



(a + l)x 

(a + 1) [(2a + l)x 2 -l 



1 

45 



a(a + l)x (2a + l)x 2 - 3 



1 



1080 
1 

18900 
1 

907200 



a(a + 1) (2a - l)(2a + l)x 4 - 6(2a - l)x 2 + 3 



(a - l)a(a + l)x (2a - l)(2a + l)x 4 - 10(2a - l)x 2 + 15 
(a - l)a(a + 1) [(2a - 3)(2a - l)(2a + l)x 6 + 



15(2a - 3)(2a - l)x 4 + 45(2a - 3)x 2 - 15 



28576800 



[a - 2)(a - l)a(a + l)x (2a - 3) (2a - l)(2a + l)x 6 + 



21(2a - 3)(2a - l)x 4 + 105(2a - 3)x 2 - 105 



So we might guess that 

4(-l) i+1 (a + l' 



c* 2i (a,x) 
c\(a, x) = 
c 2i+1 {a,x) = 



(2i)! 



i - 1 



2^1 



-i+l,a + 



x 



i = 1,2,3, 



-1) 



a 



I (a + l)x 2F1 



(2t + l)! 

Hence, we have found the following conjecture. 



-i + l,a + | — % 

3 
2 



Conjecture. The polynomials ^P" ,a,M,M (a 
the form 



(16) 



n=0 



i = 1,2,3,. 



satisfy a linear differential equation of 



M]T c;(x)y (i) (x) + (1 - x 2 )y"(x) - 2(a + l)xy'(x) + n(n + 2a + l)y{x) = 0, 



8=0 



where the coefficients {q(x)}°^ are given by (53} , ( [L5[ ) and (M). 

We remark that this conjecture would imply that for nonnegative integer values of a the 
polynomials \ P" ,a,M,M (x) \ satisfy a linear differential equation of order 2a + 4. 

I. > n=0 



Finally, we note that since 



2F1 



-i + 1, a + 



(-a-2 + »; 



i-l 



(I) 



, 1 = 1,2,3,... 



i-l 
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and 



we have 



2^1 



-i + 1, a + § — i 



1 



(-a - 1 + 

(§)* 



i 



i = 1,2,3,... 



,2^-1 

4 (-a- l)j_i (-a -2 + i)j_i 



8=1 



l(i)i(l)i4< (t-1)! 



i=0 



i!(|) i (2)i(^) < V ^ (2i)!(3) 2 , 



and 



XI c 2i+i («>-!) 

8=0 



E 



8(a + 1) (-a)t-i (-a - 1 + 
J (l)<(§)i4< 

'3 (a + ij ^ i!(2)«(§) 4 (§), ( 4 j 



"o(« + l)E 



(~«)2i 2i _ (~° ~ 1 )2i+l _ 2 2i+l 



2E 



i=0 



(2i + l)!(4) 2i ^(2t + l)!(3) 2i+1 



In the same way we find 



OO OO / -i \ 



8=1 



^ (2i)!(3) M 



and 



This implies that 



OO OO / -I \ 

E4 +1 («,i) = -2E i: a :. Wl * 2i+1 



i=0 



i=0 



(2i + l)!(3) 2 



i+l 



i=l 



i=l 
oo 



8=0 



and 



^c*(a,l) = 53<^(a,l) +E c 2i+l( a ' 1 ) 



8=1 



8=1 

OO 



8=0 



2 E ( S 5 ^ i (-2)* = 2iJi 



8=0 



i!(3)< 



-a — 1 
3 



Remark. The proof of the conjecture has been found and will be given in a forthcoming 
paper [§. 

Acknowledgement. The author wishes to thank Professors Desmond Evans, Norrie 
Everitt and Lance Littlejohn for their invitation to come to Cardiff, where the formulas of 
the conjecture were found. 

Further he wishes to thank his uncle Jan Koekoek for his valuable remarks after reading 
the first version of this paper. 



23 



References 

[1] W.N. Everitt and L.L. Littlejohn : Orthogonal polynomials and spectral theory : a 
survey. Proceedings of the Third International Symposium on Orthogonal Polynomials 
and Their Applications, Erice (Italy), 1990, to appear. 

[2] J. KOEKOEK and R. Koekoek : On a differential equation for Koornwinder's gen- 
eralized Laguerre polynomials. Proceedings of the American Mathematical Society, to 
appear. 

[3] R. Koekoek : Generalizations of the classical Laguerre polynomials and some q- 
analogues. Thesis, 1990. 

[4] R. Koekoek : A differential equation for symmetric generalized ultraspherical Jacobi 
polynomials. In preparation. 

[5] R. Koekoek and H.G. Meijer : A generalization of Laguerre polynomials. SIAM 
Journal on Mathematical Analysis, to appear. 

[6] T.H. KOORNWINDER : Orthogonal polynomials with weight function (1 — x) a (l + x)@ + 
M8{x + 1) + NS(x - 1). Canadian Mathematical Bulletin 27(2), 1984, 205-214. 

[7] A.M. Krall : Orthogonal polynomials satisfying fourth order differential equations. 
Proceedings of the Royal Society of Edinburgh A87, 1981, 271-288. 

[8] A.M. Krall and L.L. Littlejohn : On the classification of differential equations 
having orthogonal polynomial solutions II. Annali di Matematica Pura ed Applicata (4) 
149, 1987, 77-102. 

[9] H.L. Krall : Certain differential equations for Tchebycheff polynomials. Duke Mathe- 
matical Journal 4, 1938, 705-718. 

[10] H.L. Krall : On orthogonal polynomials satisfying a fourth order differential equation. 
The Pennsylvania State College Studies, No. 6, 1940. 

[11] L.L. Littlejohn : The Krall polynomials: A new class of orthogonal polynomials. 
Quaestiones Mathematicae 5, 1982, 255-265. 

[12] L.L. Littlejohn : On the classification of differential equations having orthogonal 
polynomial solutions. Annali di Matematica Pura ed Applicata (4) 93, 1984, 35-53. 



24 



